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Abstract 

This  study  investigates  the  characteristics  of  transport  phenomena  in  the  cathode  gas  diffusion  layer  (GDL)  of  a  proton  exchange  membrane 
fuel  cell  (PEMFC)  and  their  influences  on  cell  performance  by  utilizing  a  two-phase  flow  model  based  on  a  multiphase  mixture  formulation.  The 
profiles  of  liquid  water  saturation  and  oxygen  concentration  across  the  GDL  are  obtained  by  a  direct  numerical  procedure,  and  the  corresponding 
polarization  curves  are  determined  as  well.  The  results  show  that  the  liquid  water  generated  by  the  electrochemical  reaction  could  significantly 
reduce  the  effective  porosity  of  the  GDL  under  high  current  density  conditions  to  hinder  oxygen  transport  through  the  GDL.  The  removal  rate 
of  liquid  water  due  to  capillarity-induced  motion  can  be  enhanced  conspicuously  by  an  increase  in  hydrophobicity  at  first,  and  this  effect  then 
diminishes  gradually.  The  effects  of  porosity  and  thickness  of  the  GDL,  inlet  velocity  of  reactant  gas,  and  gas  channel  length  on  the  cell  performance 
are  also  examined  in  detail.  The  results  provide  insights  into  the  transport  mechanisms  in  the  GDL  and  benefit  the  design  of  PEMFCs. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is  con¬ 
sidered  as  an  alternative  power  source  for  widely  promising 
applications  in  many  power  systems.  However,  due  to  its  lower 
operating  temperature  and  the  utilization  of  external  humidifica¬ 
tion  of  the  reactant  gases,  water  management  is  quite  important 
to  avoid  the  so-called  “flooding”  phenomenon,  particularly  on 
the  cathode  side  where  the  water  vapor  may  condense  and  block 
the  pores  of  the  GDL,  and  so  reduce  the  amount  of  oxygen 
transport  to  the  catalyst  layer.  Such  a  limiting  factor  in  PEMFC 
performance  is  well-known,  and  thus  it  is  vital  to  understand  the 
two-phase  transport  characteristics  within  the  cathode  GDL  and 
their  influences  on  the  cell  performance. 

Numerous  studies  have  devoted  to  developing  mathemati¬ 
cal  models  for  the  transport  of  reactants  and  product  water  in 
a  PEMFC  [1-11].  Particularly,  since  it  involves  a  complicated 
two-phase  transport  process  in  the  GDL,  several  two-phase  flow 
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models  for  the  cathode  of  PEMFC  have  been  proposed  to  sim¬ 
ulate  the  effects  of  the  gas  and  liquid  water  hydrodynamics  on 
cell  performance.  He  et  al.  [12]  developed  a  two-phase  flow 
model  for  PEMFCs  with  interdigitated  flow  fields.  In  this  model, 
they  considered  the  transport  of  liquid  water  through  the  com¬ 
pletely  wetted  GDL  which  is  driven  by  the  shear  force  of  gas 
flow  and  capillary  force.  Berning  and  Djilali  [13]  presented  a 
3D,  multiphase,  and  multicomponent  model  for  the  anode  and 
cathode  of  a  PEMFC.  They  described  the  two-phase  flow  inside 
the  GDL  by  the  unsaturated  flow  theory  (UFT)  that  a  uniform 
gas-phase  pressure  is  assumed  within  the  GDL.  In  contrast  to 
the  traditional  UFT,  Wang  and  Cheng  [14,15]  proposed  a  multi¬ 
phase  mixture  (M2)  model  that  the  assumption  of  constant  gas 
pressure  is  released  and  it  accounts  for  the  gas  flow  counter 
to  the  capillarity-driven  liquid  flow.  This  M2  model  has  been 
widely  used  recently  to  investigate  the  multiphase  and  mul¬ 
ticomponent  transport  in  cathode  GDL  [16-21].  Wang  et  al. 
[16]  performed  a  study  for  the  two-phase  transport  of  reac¬ 
tants  and  products  in  the  cathode  of  a  PEMFC.  They  considered 
single  and  two-phase  regimes  of  water  distribution  and  classi¬ 
fied  them  by  a  threshold  current  density  corresponding  to  the 
emergence  of  liquid  water  at  the  membrane/cathode  interface. 


M.-H.  Chang  et  al.  /  Journal  of  Power  Sources  160  (2006)  268-276 


269 


Natarajan  and  Nguyen  [17]  developed  a  2D  two-phase  tran¬ 
sient  model  for  the  cathode  of  a  PEMFC  using  conventional  gas 
distributors  and  found  that  the  liquid  water  transport  is  the  slow¬ 
est  mass  transfer  phenomenon  in  the  cathode  which  is  mainly 
responsible  for  mass  transfer  restrictions.  A  two-phase  mixture 
model  was  also  reported  by  You  and  Liu  [18],  and  accordingly 
they  obtained  some  typical  two-phase  flow  distributions  in  the 
cathode  GDL  and  gas  channel.  They  also  discussed  the  main 
parameters  influencing  water  transport  across  the  membrane. 
Pasaogullari  and  Wang  [19]  studied  the  liquid  water  transport 
in  both  hydrophilic  and  hydrophobic  GDL  of  PEMLCs.  They 
made  some  simplifications  and  then  derived  a  one-dimensional 
analytical  solution  for  liquid  water  transport  across  the  GDL.  A 
detailed  exploration  for  the  effects  of  GDL  wettability  on  liquid 
water  transport  was  also  given  in  this  work.  Moreover,  they  used 
the  M2  model  to  investigate  the  effects  of  counter  gas  flow  with 
respect  to  the  liquid  water  flow  [20]  and  found  that  it  should 
be  a  new  oxygen  transport  mechanism.  They  further  considered 
the  two-phase  transport  in  PEMFCs  with  multilayer  diffusion 
media  [20,21]  and  examined  the  properties  of  the  micro-porous 
layer  including  average  pore  size,  wettability,  thickness,  and 
porosity. 

In  the  present  study,  a  two-phase  flow  model  based  on  the 
multiphase  mixture  concept  is  developed  to  investigate  the 
transport  characteristics  in  a  hydrophobic  GDL  of  PEMFC.  In 
addition  to  the  variations  of  liquid  water  saturation  and  oxy¬ 
gen  concentration  profiles  across  the  GDL,  the  effects  of  gas 
inlet  velocity,  gas  channel  length,  and  GDL  properties  includ¬ 
ing  wettability,  thickness,  and  porosity  on  the  polarization  curve 
of  a  PEMFC  were  also  examined  in  detail.  The  results  give  an 
overview  for  the  influences  of  two-phase  transport  in  GDL  on 
the  cell  performance. 

2.  Mathematical  model 


form  velocity  um,  where  L  is  the  gas  channel  length,  H\ 
the  gas  channel  thickness  (H2  —  H\ )  and  (H3  —  H2)  are  the 
thicknesses  of  gas  diffusion  layer  and  catalyst  layer,  respec¬ 
tively.  To  investigate  the  two-phase  transport  within  the  GDL, 
the  following  assumptions  are  made  to  simplify  the  model 
development: 

1 .  The  system  is  isothermal  and  under  steady  state. 

2.  The  gas  phase  is  an  ideal  gas  mixture  and  insoluble  in  the 
liquid  phase. 

3.  The  catalyst  layer  is  infinitely  thin  and  can  be  treated  as  a 
boundary  where  electrochemical  reaction  occurs. 

4.  The  membrane  is  fully  hydrated  by  liquid  water. 


According  to  the  multiphase  mixture  (M2)  model  and  above 
assumptions,  the  continuity  equation  for  the  mixture  in  the  GDL 
can  be  written  as 


dp 

sf  +  V.(pu)  =  0,  (1) 

ot 

where  s  is  the  porosity,  p  the  two-phase  mixture  density,  and  u 
the  superficial  two-phase  mixture  velocity.  The  species  conser¬ 
vation  equation  for  this  two-phase  flow  problem  [14]  is  given 
as 


8^(pCa)  +  V  •  (yapuCa) 
ot 


=  V  .  (spDaVCa)  +  V  . 


Y.mskDtivcf  -  vca) 

_  k 


2.1.  Governing  equations 

The  modeling  domain  for  the  cathode  GDL  is  shown  in 
Fig.  1.  The  humidified  air  enters  the  gas  channel  at  a  uni- 


where  C  is  the  species  mass  concentration,  D  the  diffusion  coeffi¬ 
cient,  s  the  phase  saturation,  and  the  superscript  a  and  subscript 
k  denote  the  species  and  phase,  respectively.  The  vector  j  is 
the  diffusive  mass  flux  and  y  is  the  advection  correction  factor 


Catalyst  Layer 


Gas  Diffusion 
Layer 


Gas  Channel 


Fig.  1.  Schematic  diagram  of  the  model  for  the  cathode  side  of  PEMFC. 
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defined  by 

Ya  =  ^2 - •  (3) 

J2PkSkCk 

k 

where  A.&  is  the  mobility  of  phase  k.  Since  this  study  is  dealing 
with  a  two-phase  flow  problem,  the  parameter  k  is  equal  to  2 
including  both  liquid  and  gas  phases.  Consequently,  Eq.  (2)  can 
be  reduced  to 

8^(pCa)  +  V  •  (yaPuCa) 
ot 

=  V  .  ( spDaVCa )  +  £V  .  [p^Df  VCf  +  Pg(l  -  s)Dag VCjf 

—  pDaVCa]  —  V  •  [(Ci  —  Cpjj].  (4) 

Note  that  here  s  denotes  the  liquid  water  saturation  that  repre¬ 
sents  the  volume  ratio  of  the  liquid  water  to  the  void  space  of 
the  GDL.  The  quantities  of  the  mixture  in  the  above  equation 
are  defined  as 


P  =  PlS  +  Pg(l  -s), 

(5) 

pCa  =  PisC?  +  pg(l  -  s)C%, 

(6) 

pu  =  pi ui  H-  PgUg, 

(7) 

pDa  =  pxsD^  +  pg(l  -  s)Dag, 

and  ya  becomes: 

(8) 

p(AiC“  +  AgC“) 

Ya  p\sC®  +  pg(l  -  s)C“’ 

(9) 

where  the  mobilities  y\  and  yg  are  defined  as 

,  K\/v\ 

k\  =  , 

kA/v\+kVg/vg 

(10) 

^  _  ^rg/% 

8  K\/v\  +  kXg/Vg 

(11) 

The  parameters  kv\  and  kvg  are  the  relative  permeability  of  the 
liquid  and  gas  phase  given  by  empirical  correlation  [22] : 

=  s3,  (12) 

krg  =  (l-s)\  (13) 


The  diffusive  mass  flux  of  liquid  phase  ji  can  be  given  as  [15]: 

K  A4A.fr 

jl  =  -^[VPC  +  (A  -  Pg)g],  (14) 

V 

where  K  is  the  absolute  permeability,  Pc  the  capillary  pressure, 
and  v  the  mixture  kinematic  viscosity  defined  by 


1 


v  = 


krl/n  +  krg/Vg 


(15) 


In  general,  the  gravity  term  in  Eq.  (14)  is  negligible  in  compar¬ 
ison  with  the  capillary  effect,  so  ji  can  be  expressed  by 


Kk  iA., 


■VP, 


Ci 


(16) 


and  the  capillary  pressure  between  the  gas  and  liquid  phases  can 
be  related  to  the  liquid  phase  saturation  [15]  as 

Pc  =  acos(0c)[  —  j  J(s ),  (17) 

where  cr  is  the  surface  tension  of  liquid  water,  0C  the  contact 
angle,  and  J(s )  is  the  Leverett  function  [19].  For  hydrophilic 
GDL  (0°  <0C<  90°),  J(s)  is  given  by 

./(.S’)  =  1.417(1  -s)-  2.12(1  -  s)2  +  1.263(1  -  sf,  (18) 

while  for  hydrophobic  GDL  (90°  <  0C  <  180°),  J(s )  is 

J(s)  =  1.417s  —  2.12s2  +  1.263s3.  (19) 


Substituting  Eq.  (17)  in  Eq.  (16),  we  can  obtain: 


Jl  = 


KX  iA., 


v 


CT  COS(0e)(-|:) 


£  \  !/2d  J(s) 
ds 


Vs. 


(20) 


2.2.  Water  transport 


When  the  vapor  partial  pressure  exceeds  the  corresponding 
saturation  vapor  pressure,  the  vapor  condenses  and  the  concen¬ 
tration  of  water  vapor  in  the  gas  mixture  is 


PVMU2° 

PMd" 


(21) 


Assuming  that  oxygen  and  nitrogen  are  insoluble  in  the  liquid 
phase,  the  concentration  of  liquid  water  is 

cf  2°  =  1 .  (22) 


Since  the  GDL  is  isothermal  and  under  steady  state,  if  the  GDL 
is  assumed  to  be  saturated  with  water  vapor,  the  concentration 
of  water  vapor  is  uniform  and  the  gas  phase  diffusion  vanishes. 
Accordingly,  substituting  Eqs.  (21)  and  (22)  into  Eq.  (4)  yields: 


V  •  (KH2opuCH2°)  +  V 


(23) 


where  the  advection  correction  factor  yn2o  becomes: 


Kh2o  = 


AiqH2°  +  AgCj120 
CH-° 


(24) 


Integrating  the  continuity  equation  (1)  for  the  two-phase  mixture 
over  the  GDL  [20],  we  can  obtain: 


pv  =  — 


I 

IF 


(1+2  a)MU2°  + 


(25) 


where  I  is  the  current  density,  F  the  Faraday  constant,  M  the 
molecular  weight,  and  a  the  net  water  transport  coefficient  across 
the  membrane.  A  positive  a  means  that  the  net  water  transport 
is  from  the  anode  side  to  the  cathode  side.  Combining  Eqs.  (24) 
and  (25)  we  have 


Kh2oPvCH2°  = 


+  Xg 


p^\ 

Pg  J 


x 


I 

2F 


(1  +2a)MU2°  + 


7-M°2 

4  F 


(26) 


V 
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Substituting  Eqs.  (20)  and  (26)  into  Eq.  (23)  and  then  integrating 
the  resultant  equation  over  the  GDL,  we  can  obtain  the  governing 
equation  for  the  liquid  water  saturation: 


Xg 


Pf-° 


Pg 


— —  (1  +  2a)  MHl°  +  —  M°2 
IF  4  F 


P 


h2o 


+  ,  _  ^  ,  ^co  , L y/^Vj 


Pg 


V 


\kJ 


ds 


(27) 


2.3.  Oxygen  transport 

It  is  assumed  that  oxygen  is  insoluble  in  the  liquid  phase,  so 
we  have 

cP2  =  0.  (28) 

Consequently,  Eq.  (4)  for  oxygen  becomes: 

V  •  (xo2puC°2)  =  V  •  [epg(l  -  s)D°2VC°2]  +  V  •  (C®2*), 

(29) 


liquid  water  and  the  liquid  water  saturation  is  zero  at  the  interface 
between  the  gas  channel  and  the  GDL,  we  have 

a  =  0,  aty  =  H\.  (34) 


Furthermore,  at  this  interface  the  mass  flux  of  oxygen  can  be 
written  in  the  form: 


(35) 


where  hm  is  the  convective  mass  transfer  coefficient  and  C®2  is 
the  local  average  concentration  of  oxygen.  Note  that  the  param¬ 
eter  hm  can  be  determined  from  the  Sherwood  number  defined 
by  [16]: 


(36) 


and  the  average  oxygen  concentration  decreases  along  the  chan¬ 
nel  due  to  electrochemical  reaction.  The  mass  flux  of  oxygen  is 
assumed  a  constant  along  the  channel.  Accordingly,  by  using  the 
approach  of  control  volume  integration  within  the  dotted  domain 
as  shown  in  Fig.  1,  we  can  obtain: 

■  _n  n  M°2I 

Min  Hi  pf(Cf  -  Cp?n)  =  -  L 1 .  (37) 


where  the  advection  correction  factor  yo2  is 

AgCg°2 

y°2  =  -cor- 

Combining  Eq.  (30)  with  Eq.  (25),  we  have 


(30) 


Combining  Eqs.  (35)  and  (37)  yields  the  boundary  condition  for 
oxygen  concentration  at  the  interface  y  =  H\ : 


M°2 1  /  1  Li  \ 
4 pfF  \hm  umH\) 


(38) 


Yo2pvC°2  =  XgCf 


I 

2F 


(1  +  2  a)MU2°  + 


3—M° 2 

4  F 


(31) 


Substituting  Eqs.  (20)  and  (31)  into  Eq.  (29)  and  then  integrating 
the  resultant  equation  over  the  GDL  thickness,  one  has 

— —  (1  +  2g9MH2°  +  —  M°2 
2  F  AF 

n  Kkyka  (  S\l/2dJ(s ) 

-  [c?-^acoS<(W(-) 

-  [epg(l  -  s)Dp2’cffVCp2]  =  2-M°2.  (32) 

In  above  equation,  the  effective  diffusion  coefficient  of  oxygen 
D® 2,eff  is  modified  by  Bruggeman-type  correlation  [23]  to  con¬ 
sider  the  tortuosity  factor: 


AgCg°2 


2.5.  Polarization  curve 


The  electrochemical  reaction  on  the  catalyst  layer  is 
described  by  the  Tafel  equation  in  the  form: 


/  =  (  1 


(39) 


where  7ref  is  the  exchange  current  density,  ac  the  cathode  transfer 
coefficient,  F  the  Faraday  constant,  R  the  universal  gas  constant, 
T  the  absolute  temperature,  and  rj  is  the  overpotential.  Here  the 
term  (1  —  s)  is  used  to  account  for  the  reduction  of  active  surface 
due  to  liquid  water  coverage  of  catalyst  particles.  Hence,  the 
output  cell  voltage  can  be  determined  by 


V  —  Eoc  0  Oohmi 


(40) 


Dg2’efi  =  [e(l  -  s)\TD®2.  (33) 

Eq.  (32)  is  the  governing  equation  of  oxygen  concentration.  As 
shown,  the  transport  of  oxygen  involves  the  advection,  capillary 
force,  and  diffusion  effects. 

2.4.  Boundary  conditions 

Both  Eqs.  (27)  and  (32)  are  first  order  nonlinear  ordinary 
differential  equation  and  two  boundary  conditions  are  needed 
to  obtain  the  solution.  Assuming  the  gas  channel  is  free  from 


where  Voc  is  the  open  circuit  voltage,  and  r]ohm  is  the  ohmic 
losses.  Assume  that  the  ohmic  losses  are  primary  caused  by  the 
ion  transport  through  the  membrane  and  the  membrane  is  fully 
saturated  with  liquid  water,  thus  ijohm  can  be  approximated  by 

IHm 

Oohm  =  ~  ?  (41) 

where  Hm  and  crm  are  the  thickness  and  ion  conductivity  of  the 
membrane,  respectively. 

Eqs.  (27)  and  (32)  together  with  the  boundary  conditions  (34) 
and  (38)  can  be  solved  by  a  fourth-order  Runge-Kutta  method 
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Table  1 


Physical  properties  and  parameters 


Cell  temperature,  T  (K) 

353 

Air  pressure,  Pc  (Pa) 

1.013  x  105 

Universal  gas  constant,  R  (J  mol-1  K-1) 

8.314 

Faraday  constant,  F 

96487 

GDL  porosity,  s 

0.5 

Oxygen  molecular  weight,  (kg  mol-1) 

0.032 

Water  molecular  weight,  Mh2o  (kg  mol-1) 

0.018 

Oxygen  diffusion  coefficient,  D®2  (m2  s-1) 

1.805  x  10"5 

Oxygen  mass  fraction  at  the  inlet,  C% 

0.21 

Net  water  transport  coefficient,  a 

0.5 

GDL  thickness  (H2  —  Hi)  (m) 

3.0  x  10-4 

Channel  thickness,  H 1  (m) 

1.0  X  10“3 

Liquid  water  density,  p\  (mkg-3) 

974.85 

Liquid  water  kinematic  viscosity,  v\  (m2  s-1) 

3.65  x  10-7 

GDL  permeability,  K  (m2) 

5.0  x  10-13 

Surface  tension,  a  (N  m-1) 

0.0625 

Cathode  transfer  coefficient,  ac 

1 

Inlet  velocity,  um  (ms-1) 

0.4 

Channel  length,  L  (m) 

0.05 

with  adaptive  step  size.  Some  typical  parameters  used  in  this 
work  are  given  in  Table  1 . 

3.  Results  and  discussion 

The  developed  two-phase  flow  model  for  cathode  GDL  is 
used  to  examine  the  effects  of  many  important  factors  on  cell  per¬ 
formance  including  the  GDL  contact  angle,  porosity,  thickness, 
gas  inlet  velocity,  and  channel  length.  Note  that  the  contact  angle 
0C  implies  the  wettability  of  the  GDL  and  a  larger  0C  indicates  the 
GDL  has  a  stronger  hydrophobicity.  Fig.  2  shows  the  variation 
of  profile  of  liquid  water  saturation  across  the  GDL  thickness 
with  several  assigned  values  of  0C  under  1=  14,000  A  m-2.  Obvi¬ 
ously,  the  liquid  water  saturation  reduces  significantly  for  a  more 
hydrophobic  GDL  especially  at  the  interface  between  the  GDL 
and  the  catalyst  layer  (CL).  This  effect  decays  gradually  with 


Fig.  2.  Profiles  of  liquid  water  saturation  s  across  the  GDL  for  several  assigned 
contact  angles  6C  with  /=  14,000  Am-2. 


Fig.  3.  Profiles  of  oxygen  concentration  across  the  GDL  for  several  different 
contact  angles  6C  and  single-phase  flow  with  I-  14,000  Am-2. 

increase  of  6C.  For  example,  at  the  GDL/CL  interface,  the  value 
of  s  reduces  from  0.278  to  0.126  as  0C  increases  from  91°  to 
100°.  The  percent  rate  of  change  of  liquid  water  saturation  is 
about  54.7%.  But  when  0C  increases  further  from  100°  to  110°, 
the  value  of  s  is  0. 104  at  0C  =  1 10°  and  the  corresponding  percent 
rate  of  change  of  s  is  about  17.5%  only.  Therefore,  it  should  be 
an  efficient  way  to  improve  the  flooding  phenomena  by  increas¬ 
ing  the  GDL  contact  angle,  but  this  effect  is  limited  when  0C 
is  large  enough.  The  variation  of  profile  of  oxygen  concentra¬ 
tion  is  shown  in  Fig.  3  for  several  assigned  values  of  0C  with 
1=  14,000  A  m~2.  It  is  noted  that  the  curve  for  single-phase  flow 
corresponds  to  the  state  where  there  is  no  liquid  water  within 
the  GDL.  As  seen,  oxygen  concentration  at  the  GDL/GL  inter¬ 
face  is  lower  in  the  case  0C  =  91°  because  the  higher  liquid  water 
saturation  as  shown  in  Fig.  2  blocks  the  pores  of  GDL  resulting 
in  an  increase  of  oxygen  transport  resistance.  As  0C  increases, 
the  removal  rate  of  liquid  water  rises  since  the  capillary  force  of 
GDL  becomes  greater.  Thus,  the  resistance  of  oxygen  transport 
reduces  and  the  oxygen  concentration  at  the  GDL/CL  interface 
increases  with  0C.  Such  an  effect  is  apparent  at  first  when  0C 
increases  from  91°.  As  0C  >  100°,  the  profile  of  oxygen  concen¬ 
tration  gradually  approaches  the  curve  of  the  single-phase  flow 
and  the  rise  of  oxygen  concentration  by  using  a  GDL  with  higher 
6>c  becomes  limited.  Fig.  4  illustrates  the  polarization  curves  for 
several  assigned  values  of  0C.  It  is  found  that  the  influence  of  liq¬ 
uid  water  saturation  on  cell  performance  is  not  significant  when 
the  cell  operates  under  lower  current  density  conditions.  How¬ 
ever,  at  higher  current  densities,  more  oxygen  is  consumed  and 
more  water  is  produced,  causing  the  influence  of  GDL  hydropho¬ 
bicity  to  be  pronounced.  In  the  case  0C  =  91°,  the  higher  liquid 
water  saturation  makes  the  mass  transport  losses  occur  earlier 
resulting  in  a  lower  limiting  current  density.  Although  the  limit¬ 
ing  current  density  can  be  raised  by  increasing  0C,  it  is  still  less 
than  that  in  the  case  of  single-phase  flow  due  to  the  existence 
of  liquid  water  within  the  GDL.  Furthermore,  the  increments  of 
the  limiting  current  density  between  the  cases  0C  =  100°,  110°, 
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Fig.  4.  Polarization  curves  for  several  different  contact  angles  0C  and  single¬ 
phase  flow. 

and  120°  are  quite  small,  which  indicates  the  rate  of  oxygen 
transport  under  such  a  high  current  density  is  almost  invariant 
with  0C  when  0C  is  sufficiently  large. 

The  effects  of  GDL  porosity  on  the  profile  of  liquid  water 
saturation  are  shown  in  Fig.  5  for  three  typical  values  of  e.  Both 
cases  0C  =  91°  and  110°  are  considered  with  I-  13,000  Am-2. 
For  0C  =  9 1  ° ,  it  is  found  that  the  liquid  water  saturation  decreases 
with  an  increase  of  GDL  porosity,  which  indicates  that  a  larger 
porosity  would  enhance  the  removal  rate  of  liquid  water  in  the 
GDL.  But  this  effect  diminishes  for  a  more  hydrophobic  GDL  as 
shown  in  the  case  0C  =  1 10°.  Therefore,  the  porosity  may  influ¬ 
ence  the  liquid  water  saturation  profile  significantly  only  when 
the  GDL  has  a  lower  hydrophobicity. 


Fig.  5.  Profiles  of  liquid  water  saturation  s  across  the  GDL  for  both  cases  of  0C 
with  three  assigned  GDL  porosities  under  1=  13,000  Am-2. 


Fig.  6.  Profiles  of  oxygen  concentration  across  the  GDL  for  both  cases  of  6C  and 
single-phase  flow  with  three  assigned  GDL  porosities  under  1=  13,000  Am-2. 

In  contrast  to  the  small  influence  on  liquid  water  saturation 
profile,  the  GDL  porosity  plays  an  important  role  in  the  varia¬ 
tion  of  oxygen  concentration  profile  as  demonstrated  in  Fig.  6. 
A  larger  porosity  significantly  reduces  the  resistance  of  oxy¬ 
gen  transport  through  the  GDL  and  greatly  raises  the  oxygen 
concentration  on  the  GDL/CL  interface  due  to  the  increase  of 
effective  diffusion  coefficient  as  indicated  in  Eq.  (33).  For  exam¬ 
ple,  the  oxygen  concentration  adjacent  to  the  catalyst  layer  in 
the  case  s  =  0.3  with  0C  =  91°  is  almost  consumed  thoroughly  by 
the  electrochemical  reaction.  The  insufficient  supply  of  oxy¬ 
gen  causes  severe  concentration  losses  and  restricts  the  cell 
performance.  Such  phenomena  can  be  improved  if  the  GDL 
material  is  replaced  by  a  higher  porosity  one  as  shown  in  the 
cases  £  =  0.4  and  0.5.  Similar  results  also  can  be  observed  in 
the  cases  0C  =  110°  and  the  single-phase  flow  model.  It  also  can 
be  predicted  from  these  results  that  the  limiting  current  density 
can  be  enhanced  by  a  larger  porosity  for  a  given  GDL  contact 
angle  Qc.  Fig.  7  shows  the  corresponding  polarization  curves  for 
these  three  typical  porosities  s  =  0.3,  0.4,  and  0.5,  and  each  case 
contains  three  typical  cases  0C  =  91°,  110°,  and  the  model  of 
single-phase  flow.  Apparently,  the  effects  of  two-phase  flow  are 
important  when  the  fuel  cell  operates  under  high  current  den¬ 
sity  conditions  and  the  cell  performance  is  better  with  higher 
porosity. 

The  thickness  of  cathode  GDL  is  also  an  important  factor 
to  affect  the  cell  performance.  Fig.  8  illustrates  the  profiles  of 
liquid  water  saturation  for  three  different  thicknesses  200,  300, 
and  400  jam.  Both  typical  contact  angles  0C  =  91°  and  110°  are 
considered  with  current  density  1=  14,000 Am-2.  The  results 
show  that  if  the  GDL  is  less  hydrophobic  (i.e.  0C  =  9 1  °),  the  lower 
removal  rate  of  liquid  water  causes  a  high  liquid  water  saturation 
on  the  GDL/CL  interface.  Moreover,  a  thicker  GDL  leads  to  a 
higher  resistance  for  liquid  water  transport.  As  a  result,  the  liquid 
water  saturation  rises  rapidly  with  increasing  GDL  thickness. 
However,  if  the  GDL  is  more  hydrophobic  (i.e.  0C  =  110°),  the 
capillarity-induced  flow  dominates  the  liquid  water  transport 
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Fig.  7.  Variation  of  polarization  curves  with  contact  angle  0C  for  three  assigned 
GDL  porosities:  (a)  s  =  0.3,  (b)  £  =  0.4,  and  (c)  e  =  0.5. 

and  the  removal  rate  is  increased.  Therefore,  the  liquid  water 
saturation  on  the  GDL/CL  interface  reduces  and  the  variation  of 
GDL  thickness  only  produces  a  slight  influence  on  the  profile  of 
liquid  water  saturation. 

The  effects  of  GDL  thickness  on  the  profile  of  oxygen  concen¬ 
tration  are  demonstrated  in  Fig.  9  for  the  same  cases  considered 
in  Fig.  8  and  the  case  of  single-phase  flow.  As  shown  in  this 
figure,  the  oxygen  transport  is  profoundly  affected  by  the  GDL 


Fig.  8.  Profiles  of  liquid  water  saturation  s  across  the  GDL  for  both  cases  of  6C 
with  three  assigned  GDL  thicknesses  under  I-  14,000  Am-2. 


Fig.  9.  Profiles  of  oxygen  concentration  across  the  GDL  for  both  cases  of  6C  and 
single-phase  flow  with  three  assigned  GDL  thicknesses  under  I-  14,000  A  m-2. 

thickness.  In  the  case  0C  =  91°,  the  oxygen  concentration  near  the 
catalyst  layer  decreases  rapidly  with  an  increase  of  GDL  thick¬ 
ness.  Consequently,  the  concentration  losses  due  to  insufficient 
mass  transport  of  oxygen  take  place  in  advance  for  a  thicker 
GDL.  Similar  results  also  can  be  observed  in  the  cases  0C  =  110° 
and  the  single-phase  flow  model.  An  oxygen  concentration  dif¬ 
ference  is  noticed  on  the  GDL/CL  interface  between  the  case 
0c  =  9l°  and  the  single-phase  flow  model,  and  this  difference 
gradually  grows  with  GDL  thickness.  This  result  indicates  that 
the  influence  of  two-phase  flow  is  less  important  for  a  thinner 
GDL,  while  it  becomes  a  critical  factor  for  a  thicker  GDL  in 
fuel  cell  performance.  The  corresponding  polarization  curves 
for  these  three  typical  GDL  thicknesses  are  shown  in  Fig.  10. 
Apparently,  the  cell  performance  is  lower  for  a  thicker  GDL 
since  the  limiting  current  density  decreases  with  GDL  thick¬ 
ness. 

The  influences  of  gas  inlet  velocity  and  gas  channel  length 
also  can  be  examined  via  the  developed  theoretical  model. 
Fig.  11  demonstrates  the  variation  of  polarization  curves  for 
three  assigned  inlet  velocities.  In  the  case  %  =  0.2ms_1,  the 
volumetric  gas  flow  rate  is  too  low  to  provide  sufficient  oxygen 
through  the  GDL  for  electrochemical  reaction.  Therefore,  the 
concentration  losses  occur  early  and  make  the  lower  limiting 
current  density.  In  addition,  the  generation  of  liquid  water  is  too 
slow  to  produce  a  significant  effect  on  blocking  the  GDL.  Hence, 
the  liquid  phase  plays  a  less  important  role  when  the  gas  inlet 
velocity  is  small,  which  can  be  interpreted  by  the  little  differ¬ 
ence  between  the  limiting  current  densities  of  the  curves  for  the 
case  0C  =  91°  and  the  single-phase  flow  model.  As  um  increases, 
more  oxygen  can  be  transported  to  the  catalyst  layer  for  elec¬ 
trochemical  reaction.  As  a  result,  the  cell  performance  enhances 
and  the  limiting  current  density  rises  obviously  as  shown  in  the 
case  um  =  0.4  ms-1.  Furthermore,  the  influence  of  liquid  water 
grows  gradually  when  the  cell  operates  under  high  current  den¬ 
sities  especially  for  a  less  hydrophobic  GDL.  One  can  observe 
that  the  difference  of  limiting  current  density  between  the  case 
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Fig.  10.  Variation  of  polarization  curves  with  contact  angle  6C  for  three 
assigned  GDL  thicknesses:  (a)  H2~H\  =  200  |xm,  (b)  H2—H1-  300  pan,  and 
(c)  400  [xm. 

0C  =  91°  and  the  single-phase  model  becomes  distinct,  and  such 
a  difference  enlarges  further  in  the  case  U[n=0.6  ms-1. 

Fig.  12  shows  the  variation  of  polarization  curves  for  three 
different  gas  channel  lengths.  For  the  shorter  channel  length  case 
with  L-  5  cm,  the  fuel  cell  exhibits  a  better  performance  due  to 
the  higher  average  oxygen  concentration  along  the  channel  as 
well  as  the  electrochemical  reaction  rate.  The  higher  limiting 


Fig.  11.  Variation  of  polarization  curves  for  three  assigned  inlet  velocities  with 
both  contact  angles  6C  and  single-phase  flow. 


Fig.  12.  Variation  of  polarization  curves  for  three  assigned  gas  channel  lengths 
with  both  contact  angles  6C  and  single-phase  flow. 

current  density  simultaneously  accompanies  a  faster  water  gen¬ 
eration  rate  resulting  in  higher  liquid  water  saturation  in  the 
GDL.  Accordingly,  the  effective  porosity  of  the  GDL  reduces  as 
well  as  the  oxygen  transport,  so  the  limiting  current  densities  of 
both  cases  0C  =  91°  and  110°  are  distinctly  less  than  that  of  the 
single-phase  flow  model.  As  the  channel  length  increases,  the 
mean  oxygen  concentration  along  the  channel  descends  quickly 
which  causes  a  great  reduction  of  limiting  current  density  as 
shown  in  the  cases  L  =  1  and  14  cm.  The  effect  of  liquid  water  sat¬ 
uration  also  decreases  gradually.  For  L-  14  cm,  the  three  curves 
almost  overlap,  which  indicates  the  liquid  water  saturation  is 
quite  small  and  the  results  are  equivalent  to  those  of  the  single¬ 
phase  flow  model. 

4.  Conclusions 

A  two-phase  flow  analysis  has  been  implemented  to  exam¬ 
ine  the  characteristics  of  transport  phenomena  in  the  cathode 
GDL  of  a  PEMFC  and  their  influences  on  cell  performance.  The 
developed  theoretical  model  provides  a  simple  way  to  estimate 
the  effects  of  the  gas  inlet  velocity  and  gas  channel  length.  It 
was  known  that  a  hydrophobic  GDL  possesses  a  higher  removal 
rate  for  liquid  water  to  prevent  the  flooding  phenomenon.  There¬ 
fore,  this  study  considered  the  hydrophobic  GDL  with  0C  >  90° 
only  and  found  that  the  liquid  water  saturation  indeed  is  reduced 
significantly  under  high  current  density  conditions  when  0C 
increase  from  91°  and  the  oxygen  transport  is  enhanced  due 
to  the  improvement  of  effective  porosity  of  the  GDL.  However, 
this  effect  decays  gradually  with  0C  and  is  limited  if  0C  is  large 
enough.  The  profiles  of  liquid  water  saturation  and  oxygen  con¬ 
centration  are  sensitive  to  the  porosity  and  thickness  of  GDL.  It 
is  found  that  a  smaller  porosity  and  a  thicker  GDL  produce  the 
same  effect  to  raise  the  resistance  of  oxygen  transport  and  lower 
the  cell  performance. 

The  gas  inlet  velocity  and  the  channel  length  are  both  impor¬ 
tant  factors  that  affect  the  oxygen  concentration  along  the  chan- 
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nel  as  well  as  in  the  GDL.  The  results  show  that  the  concentration 
losses  are  serious  for  a  smaller  inlet  velocity  or  a  longer  chan¬ 
nel  length.  The  generation  of  liquid  water  is  not  significant  in 
these  cases  due  to  the  lower  limiting  current  densities.  How¬ 
ever,  an  increase  of  inlet  velocity  or  a  shorter  channel  length  can 
enhance  the  cell  performance  obviously,  and  the  two-phase  flow 
in  GDL  exhibits  an  important  influence  on  the  variation  of  lim¬ 
iting  current  density.  These  results  are  beneficial  to  the  design 
of  a  PEMFC  with  regard  to  water  management  problems. 
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